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ABSTRACT 

Charged partiele motion in the guiding center approximation is analyzed 
for models of the Jovian and Saturnian magnetosph^sric magnetic fields based on 
Voyager magnetometer observations. Field lines (jire traced and shown to 
exhibit the previously recognized [Connerney et *91,, 1981a, b] /^istention which 
arises from arimuthally circulating magnetospheric currents. The spatial 
dependencies of the guiding center bounce period and azimuthal drift rate are 
investigated for the model fields, the bounce period may be shorter or longer 
by a factor of 1-3 than in the field of the planetary dipole alone depending 
on whether a particle mirrors close to the magnetic equator and experiences 
predominantly an enhanced mirror force or at high latitude and is affected 
principally by the extended length of the field line, Non-dipolar effects in 
the gradient-curvature drift rate are most important at the equator and affect 
particles with all mirror latitudes. The effect is a factor of 10-15 fof’ 
Jupiter with its strong magnetodisc current and 1-2 for Saturn with its 'more 
moderate ring current. Limits of idiabaticity, where particle gyroradii 
become comparable with magnetic rcale lengths, are discussed and are shown to 
occur at quite modest kinetic energies for protons and heavier ions. 

INTRODUCTION 

Owing to the rapid rv>tation rates of Jupiter and Saturn and the large 
radial :sxtenta of their magnetospheres, diamagnetic currents carried by 
magnetospheric plasma rotating azimuthally around the planets contribute 
significantly to the magnetospheric magnetic fields. This effect is 
especially important for Jupiter with its stronger internal field and hence 
larger magnetosphere. It is an effect imbedded in pre-encounter ideas [Hines, 
1964; Gledhill, 1967; Piddington, 1969; loannidis and Brice, 19713, confirmed 
„ by the first Pioneer magnetometer measurements at Jupiter in 1973 [Smith et 
al., 1974] and Saturn in 1979 [Smith et al., 1980], and incorporated 
» immediately into post-encounter Jovian models [Barish and Smith, 1975; Beard 
and Jackson, 1976; Gleeson and Ax ford, 1976; Goertz et al., 1976]. The 
Voyager 1 and 2 missions in 1979 and 1981 have expanded enormously our 
knowledge of the Jovian and Saturnian magnetospheres. 
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In this paper we use planetary magnetic field and magnetodisc ring 
current models developed to fit Voyager vector magnetometer observations [Ness 
et al,, 1979a, b; 1981] to compute guiding center bounce periods and 
gradient-curvature drift speeds in the Jovian ^nd Saturnian magnetospheres. 

We find that magnetic field components arising from azimuthal plasma currents 
have a substantial effect on drifts, which for Jupiter are a factor of 10 or 
more larger than in a simple dipole approximation. Such drift corrections 
must be taken into account in calculating anisotropies of >100 keV electrons 
and ions in the middle magnetospheres of both planets. They also appear to be 
significant for calculating absorption lifetimes by the moons Ganymede 
(Jupiter) and Rhea (Saturn) and they should be dramatically important at 
Callisto (Jupiter), although there is no reported particle absorption th^re. 

t 

Our nterest is with magnetically trapped ions and electrons, adiabatic 
in the usual sense that the ratio c = r^/L of their gyroradii r^ to the scale 
length L of the magnetic field is small and energetic in the sense that ExB 
drift speeds u^ are of 9(e) compared with thermal speeds v. We are thus 
concerned with the situation where ExB, gradient, and curvature drifts are of 
the same order in c. 

MAGNETIC MODELS 

We use topologically similar current sheet magnetic models which have 

been developed to fit Voyager observations at Jupiter CConnerney et al., 

1981a] and Saturn CConnerney et al., 1921^]. In each instance the magnetic 

field consists of a contribution arising from processes internal to the planet 

plus the field of a current flowing azimuthally in a planet-centered 

cylindrical coordinate system with its z-axis aligned with the planetary 

dipole. In such coordinates the modv®-l magnetic fields are both azimuthally 

symmetric and north-south symmetric (i.e., with respect to the z = 0 plane). 

Magnetic field lines lie entirely in meridional planes (They have only B , B 

P 2 

components.), the observed B^ sweepback being an effect soon to be incorpora- 

¥ 

ted into such models CConnerney et al., 1982], The plasma current flows in a 
clockwise direction, when viewing along the z-axis, which is also the 
direction of corotation when the dipole and planetary rotation axes are 
aligned, as they essentially are in the case of Saturn. 
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The z' integral can also be done. The resultant expression is lengthy, so we 
adopt the following convenient notation 
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Because of the z-synunetry we limit ourselves to z. > 0 and obtain 
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z < D 


z > D 


(-1) (z I z“| ) indicates a term, identical to the entire preceding 
expression with z* replaced by j z“| (so that S^'*’ and are also replaced by 
Sg“ and s^”) and having an algebraic sign determined by the sign of z-D. The 
remaining integral is then carried out numerically. 
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At th<! radi/A distances of Interest, only the dipole term from the 
internal field f;'a important. Adding it to the current sheet contribution, we 
obtain for the field components 
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(5b) 


Since the only appearance of z in Eq. (3) is in the limits of integration, 3Ap/ 
3z is most readily obtained by differentiating this form 
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On the other hand, no such compact representation of 3A^/3p is apparent. We 
have therefore evaluated this derivative numerically. 

Shown in Fig. 1a as solid curves are Jovian magnetic field lines obtained 
by integrating the equation 


dp B 

p 



beginning at equatorial (z=0) points displaced at 5 Rj intervals from the 
center of the planet. Note the different scales on the vertical and 
horizontal axes. The 7 dotted curves are dipole field lines which intersect 
Jupiter's surface at the same latitudes as the 7 distended field lines. In 
the magnetodisc (p^ > 15 R^) the remnant magnetic field at the equator is 
weak, so that much of this equatorial region maps to a small range of magnetic 
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latitude around 73* at Jupiter’s surface. In identifying accurately the 
intersection latitudes we have used the fact that for any meridional field pA^ 
is constant along field lines. Finally for later use in specifying mirror 
points of trapped particles, there appears in Fig. la a grid of constant 
magnetic latitude lines# 

Figure lb is the corresponding situation for Saturn’s magnetosphere. 

Field lines are once again distended but to a far lesser degree than for 
Jupiter's magnetodisc. 

Although our methods for calculating the magnetic field differ, the 
results in Figs. 1 agree with those of Connerney et al, [193 la, 1981bJ. 

GRADIENT AND CURVATURE DRIFTS 

We begin with the well known, relativistically correct expression 
[Northrop, 1963] 
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for the 0(e) component of the guiding center drift velocity due to spatial 

A 

inhoraogeneity of the magnetic field. Here e- = B/IBI, v„ = v»e , v, = 

Cv*^-v„^)' , ae,/3s s e • V e, is the derivative of e, with respect to 

' 22“1 /2 ' 

distance s along field lines, y = (1 - v /c ) is the Lorentz factor, and 

and c, are the rest mass and (signed) charge of the drifting particle. 

From conservation of magnetic ^icrgnt M 


we obtain 
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Furthermore it is shown readily that 
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where the i subscript denotes the vector component perpendicular to fl. (I 
totally perpendicular to p in the models used in this paper.) Using Eqs. 
(7-9) in Eq. (6) we write 
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Note finally that in terms of the field strength B^^ at the mirror points of 
the bouncing (as well as drifting) guiding centers 
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Thus the energy dependence is a strictly multiplicative factor, while spatial 
and pitch angle dependencies, the latter entering through are in general 
intertwined in a complicated way. 

Like the current density J (cf. Eq. 1), BxvB is totally azimuthal 
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Here we have separated terras which depend on pitch angle through B,jj from those 
which do not. 

A useful measure of the azimuthal drift rate is the angular velocity « s 
Rj_/p, where u la a signed quantity, positive in the direction of increasing 
w varies on the rapid time scale of the guiding center bounce motion boch 
because of its dependence on p, which varies with position along field lines, 
and the dependence of on B, its derivatives, and X along the bounce path. 

This rapid time variation can be removed by averaging w over the bounce. We 
denote this bounce average by < « > 


. 1 ds , 1 ds R, 

< u > = — # — u s — ^ (17) 



Here the v,,"^ factor accounts for the varying amounts of time spent at 
different points s along the field line, the integral extends over a complete 
bounce path to and fro between each mirror point, and the bounce period Tg is 
the usual 


ds 



Because of the north^south symmetry assumed here, the integrals in £qs. 
(17) and (18) each reduce to 4 times the integral from the equator, which we 
take to be the s = 0 surface, to the northern mirror point, denoted by s^^j. 
Furthermore by Eqs. ( 8 ) and (11) we can relate v„(s) to B(s), B^^, and the 
magnetic moment H which is a constant of the bounce motion 
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where, by definition, 

For comparison of our result here with dipole results CHamlin et al., 
1961; taw, 1961; Thomsen and Van Allen, 198U we define 
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is th« equatorial crossing distance in planetary radii or Rg of the dipoie 
field line which would pass through the mirror point, whose coordinates are 
^m* easily shown that for any dipole field H# es defined oy Eq. 

(21), is independent of L and depends only on mirror latitude 1^^^, Thus for a 
dipole field H is independent of position along the lines of constant magnetic 
latitude shown in Figs. 1. 

Such is not the case when the magnetic fields arising from plasma 
currents are added. We have evaluated H(Pq, Xjjj) by numerically integrating 
Eq. (21) using the model field for particle motion along field lines 
intersecting the equrtor at different distances and also for different 
values of mirror latitude Given an observational position and a particle 
detector look direction, can be estimated from Figs, 1 ana can be 
determined by magnetic moment conservation. The latter requires magnetic 
intensity information which is available in the publications of Conni-rney et 
al. Cl98ia,b3, 

The H values in Figs. 2 corresponding to p^ s 0 are the L-independent 
dipole values, since as p^ 0 the magnetic model is dominated by its dipole 
component. We thus see that for particles mirroring at x^^^ > 20“ , H and hence 
tg are larger than they would be in a dipole field for a particle mirroring at 
the same The reason is that the model field line is longer than the 

corresponding dipole line because of the stretching out of magnetic field 
lines by the plasma current. This effect is much more important for Jupiter 
than it is for Saturn, as one would expect by visual inspection of Figs. 1. 

In Fig. 2b the tapering off of the slope of the 20° contours near p^ ~ '15 
^2 Is 3 result of the fact that the distention of field lines abates and they 
become more dipole-like at radial distances beyond the current sheet. 

At smaller X^, H generally has a smaller value than it does for a dipole. 
Here the difference in length between model and dipole field lines becomes 
secondary to the difference in mirror forces between them. Recall [Roederer, 
1970] that a particle mirroring near a point s^ of minimum intensity along a 
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magnt&lo fitld line U.e., a point where )B/3s s o) execute3 simple harmonic 
motion described by the equation 
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For the ranges of in Figs, 2 a single intensity minimum occurs on each 
model field line at the z s o magnetic equator with the exception of a small 
region around p^ = 15 near the outer terminus of Saturn's ring current. We 
shall return shortly to a brief discussion of this anomalous region. For any 
meridional field with continuous components and having north-south symmetry 
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with the right hand side, of course, evaluated at z = 0, The curves marked 
i = o“ in Figs, 2 have been constructed using Eqs, (25) and (26), The p = 

™ 1/p 

0, X s 0 dipole value of H = t/(18) * follows directly, The diminution from 

^ 2 
this dipole value for > 0 is due mostly to a large (aB^/az) term in Eq* 

(26) which reflects the 'hairpin* character of field lines near the equator, 

especially in the p^ > 15 Rj region of Jupiter's magnetosphere. 
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An Interesting aberration occurs in the k'* 15 Rg near equatorial region 
of Saturn’s magnetosphere. Over a small radial range of perhaps .5 Rg in 
extent, the point of minimum B bifurcates from the z = 0 plane into the 
northern and southern magnetic hemispheres. As this region is approached 
along the equator, B” diminishes so that an equatorial particle executes 
slower bounce motion which is manifested by an increasing Hg. When the point 
of minimum B leaves the equator, B" changes sign, z s 0 is no longer an 
oscillation point, and no longer has a real value. It is apparent from 
Fig. 2b that the field minimum departs less than 10*^ in latitude from the 
equator, for particles oscillating about the equator can mirror at = lO®. 
Such a non-equatorial trapping geometry was suggested for Jupiter by Barish 
and Smith [19753. At p s 16 R,, outside the Saturnian ring current, is 
again finite for all X^jj, indicating that the point of minimum B is once more 
at z = 0. 


With regard to drifts, let us again establish contact with previous 
dipolo work [Hamlin et al,, 1961; Lew, 1961; Thomsen and Van Allen, 19813 6y 
defining (cf. Eqs. (16) and (20)) 
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Here B is the magnitude of the dipole field at the planetary equator (4 G 
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for Jupiter; .209 G for Saturn) and L is as previously defined in Eq. (23). 
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with F/G independent of particle species and energy. In general F/G is a 
function of and p^, but in the case of a pure dipole field its (s L) 
dependence vanishes, and, like H previously, F/G becomes a function of g^ 
alone. 

Figures 3a and 3h show F/G for the Jovian and Saturnian models 
respectively in the same format used in Figs. 2. The L-independent dipole 
values /are again the limiting forms as p^ + 0 and the raagnetospheric current 
becomes a negligible factor. 

For both magnetospheres the most significant deviation from the dipole 
results is for equatorially mirroring = o particles). For such particles 
Ed* (28) becomes 
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which upon substitution from Eq. (16), noting that B^(ZsO) = 0, B^ = B, and 

s LR , « becomes 

0 J,s 
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(The dipole case where F/G(X^sO) = 1 follows immediately from Eq. 29*) The 
s 0® curves in Figs. 3 are plotted using Eq. 29. 


For Jupiter the deviation from the dipole result is an order of magnitude 

_2 

or more near p^ = 25 Rj* The enhancement is a consequence of the B factor 
in Eq. (29) and the very weak field at the center of the current sheet. In 


15 




// 


fact the magnetic field ia so weak and changes so rapidly due to the strong 
3Bp/a« that one should worry about the smallness of c and the validity of the 
adiabatic approximation. Knowing that one can go to low enough particle 
energy that c is small, we shall proceed here and return in the next section 
to quantify this limitation. 

We have evaluated F/G numerically and found for Jupiter very little 
dependence on at given p over the range 10“ < < 6o“ . Hence we have 

drawn in Fig. 3a only two of the closely nested curves in this interval. 
The reason is that the numerator integral in £q. 28 
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picks up a large (almost its entire) contribution from the 2=0 region where 
B is small: a particle may mirror far from the equator but the major 
contribution to its drift motion occurs in that brief interval during which it 
passes through the equator. The denominator integral 


. m , 

/ ds 
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on the other hand has contributions from all s values, with a sizable amount 
coming from the region about the mirror point. Thus LR,H, increases with X_ 

at fixed p as the particle traversal distance increases. However, in a 

® 2 
nearly compensating fashion L decreases with X^ (at fixed p^) so that L Hj is 

approximately constant. 

Note finally that (F/G) , reverses sign for X = 0 particles in the region 

iJ 01 

Pq > 30 Rj: the sense of the gradient drift changes from generally pro§rade 
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(With respect to the direction of planetary rotation'^ to retrograde, or vice 
versa, depending on the sign of the electrical charge on a particle. This 
reversal .results from a shift in the algebraic sign of 3B^/3p(z s o) . It is 
an effect which is confined to the equator and is no longer seen in (F/G)j (X^j 
= 10**, . 0 ^) : other effects taking place closely off z s 0 and which are also 
important because of weak B scon dominate. 

For Saturn the magnetic field at the z=0 center of the ring current is, 
relative to the field at higher latitudes, nowhere near as weak as at Jupiter. 
Hence equatorial drifts, while important, are not of sole significance. 

Rather than being an order of magnitude larger than dipole, (F/G)g is more 
usually a factor of 2 greater. What is evident in Fig. 3b and is not apparent 
in Fig. 3a are discontinuitias in F/G at the inner and outer edges of the ring 
current. These discontinuities result, cf. Eq. (29) i from discontinuities in 
BB^/ap as I(p). Gf. Ecj. (15), turns on and off at p^ s 8.5 a<id 15.5 Rg 
respectively. A discontinuity is also present at p^ s 5 Rj in Fig. 3a, but 
its evidence is masked by the dominance of the continuous 3Bg/3p arising from 
the Jovian dipole. The discontinuity is evident to a lesser extent at Saturn 
for Xjjj = 10", bub it clearly has disappeared by the time X^^^ = 30“ is reached. 

Were the current distribution I(p) continuous, these discontinuities 
would undoubtedly disappear. However, we expect the qualitative character of 
(F/G)g would be preserved, so long as significant 3l/3p is concentrated near 
8.5 Rg and 15.5 Rg. 

t 

Note finally the gap in the X„ s o“ curve for 15 Re < ^ 15.5 R-. 

Although Eq. (28) yields a real value for 0'/G)g in this interval, we know 
from our bounce time discussion that z s o is no longer a stable oscillation 
point: the canceling numerator and denominator integrals in Eq. (28) for z=0 
particles are in this instance both imaginary. 

DISCUSSION 

Weak B(zsO) means that particles have large gyro-radii at the equator. 
Large 3B^/3z (z=0) means a short scale length L for the magnetic field. The 
adiabatic approximation and our results are hence suspect, particularly for 
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Jupiter. As a measure of the range of validity of the adiabatic approximation 
we have determined the value of particle kinetic energy K for which c s 1 at z 
s 0. This value depends on 1^^^, for particles mirroring high off the equator 
pass through it with shallow pitch angles and consequent small gyro-radii 
compared with particles of the same energy which are mirroring at the equator. 

If we take L to be B(z=0)/3B /8z(zs0) it is straightforward to show that 

^ 2 
the c s 1 condition leads to the following equation for K = 



3xlo“^ B^(z=0) 1 

(1 + C 

m c^ 3B /az(zso) sina^ 
Op e 


32 j 1/2 


- 1 


(30) 


Here Z is the ionization state of the particle (+1 for protons, -1 for 
^ 2 

electrons etc.) and m^c its rest energy in MeV. particles with kinetic 

energy lower than that given by Eq, (30) have gyro-radii smaller than L, and 

hence this K is the upper energy limit for adiabaticity. In Fig. 4a we have 

plotted the characteristic factor = 3xl0"^ B^(z=0)/aB /az(z=0) for that 

c p 

portion of the Jovian magnetosphere where adiabaticity is suspect. 


Note that this quantity is always at least an order of magnitude larger 
2 

than me = .511 MeV for electrons, so that K » K /sine is at least 5 MeV in 
o e c e 

magnitude. We thus conclude that electrons up to approximately 5 MeV are 
adiabatic and have been treated validly in this paper. The limit is in fact 
considerably higher for electrons with larger because of the (sina^)”^ 
factor in Eq. (30). 

The case is, however, different for protons and heavier ions. Shown also 
in Fig. (4a) are solutions to Eq. (29) for protons with 3 different values of 
sin0^. Values of sinS^ = .05 and ,01 correspond roughly to particles 
mirroring at X^^^ s 30® and x^ = 60® in this p^ region of the Jovian 
magnetosphere. Note that the adiabaticity condition is quite severe for 
equatorial protons, with K (sind - 1) being as small as 15 keV. However, for 
protons mirroring at higher latitudes Kp is greater than 600 keV. 
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For htavi«r Ions and for all but the saiallest sih9 , we can expand the 

2 1/2 2 * 

square root In Eq. (30) according to (l+c*")' ^ » 1 + c /2 and obtain 


2 2 
2m„c^sin 
o e 



For s>* with sine^ s .05, is 167 keV at = 29 Rj, and thus a sizable 
portion of the sulfur ions in Jupiter's magnetosphere may be non-adiabatic. 


Figure i|b shows the corresponding situation for Saturn. Note that 
because Saturn's field lines are far less distended than Jupiter's, 

3B^/3z(z=0) is much smaller, and the maximum energy for which particles are 
adiabatic correspondingly higher. All electrons below 300 MeV are at least 
marginally adiabatic. The 3 proton curves correspond once again to particles 
mirroring at the equator ( sinews O and at approximately 30'“(sine^ = .4) and 
60®(sin9 = .05) magnetic latitude near p = 12 R„. 

In Figures 5 we plot for each magnetosphere the value of K for which <u> 
is equal in magnitude to the planetary angular rotation frequency 
S2j 2(2it/Pj a 10 hrs., Sir/Qg » 10.7 hrs.). Recall that <w> depends on the sign 
of particle charge and is retrograde with respect to planetary rotation for 
electrons and prograde for protons, expect for equatorially mirroring (x^zO®) 
particles in the p_ > 30 R , 14 Ro < < 15 R® regions where 3B,/3p 

switches sign (cf. Eq. 29) and the reverse is true for both species. If we 
neglect the nonalignment of Jupiter's dipole and rotation axes, the values in 
Fig. 5 are thus in general those at which an electron has no drift motion in a 
non-rotating coordinate system: <w> just balances the ^xB motion which drives 
particles to co-rotation. Corresponding values of have not been plotted 
but are approximately 1/2 to 3/4 as large as K^, depending on the value of p^. 

With the exception of those equatorial regions where 3B /3p is small so 
that compensatingly high particle kinetic energy is needed to produce 
significant gradient-curvature drifts, the values of in Figs. 5 are well 
within the range of validity of adiabatic theory. On the other hand, the 
situation is different for protons; values of Kp scaled from figs. 5 are of 
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the order of, and In the case of equatorial protons less than, the limits 
displayed in Figs. 4. 


CONCLUSIONS 

The adiabatic motions of charged particles in the distended Jovian and 
Saturnian magnetospheres are qualitatively similar to those in a dipole field 
but are quantitatively quite different. Guiding 'center bounce times may be 
shorter or longer depending on a particle's mirror latitude; particles 
mirroring near the magnetic equator experience a stronger-than-dipole mirror 
force and hence bounce more rapidly than in a dipole field; on the other hand, 
particles mirroring far from the equator bounce along field lines whose length 
is significantly enhanced by the distention, and hence such particles have 
bounce times longer than if they mirrored at the same spatial location in a 
dipole field. For Jupiter's magnetosphere non-oiplar effects on the bounce 
time are important for radial distances ranging from 10 Rj to 35 Rj* the 
largest distance we have considered, and they vary by a factor of roughly 3 
from dipole results. For Saturn the ring current affects the bounce over the 
6 Hg-lb Rg region, but owing to its weakness relative to Jupiter's magnetodisc 
current the dipole correction factor is no more than 1.2 - 2. 

The bounce-averaged angular drift rate is significantly affected by the 
field distention. This is so even for particles mirroring at high latitudes 
where the distention is small: much of their azimuthal drift takes place, 
especially in the case of Jupiter, during passage through the magnetic 
equator. In the case of Jupiter the angular gradient-curvature drift is as 
much as 10-15 times the dipole value and is of most importance between 20 Rj 
and 35 Rj. For Saturn the correction is more nearly a factor of 2 and is of 
most importance near the inner and outer edges of the ring current at 8.5 Rg 
and 15.5 Rg. 

The direction of the gradient-curvature drifts is generally that of 
co-rotation for protons and in the opposite sense for electrons — the same as 
it is for the planetary dipole. In the midcle magnetosphere of each planet 
the magnitude of the gradient-curvature drift in the model field is 
approximately equal in magnitude to the co-rotation speed for electrons and 
protons in the 100 keV - 1 MeV range. 
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Corrections, which we have evaluated are thus important in assessing the 
anisotropies of lOO keV and higher energy particles and should also be 
accounted for in calculating losses to moons orbiting through the middle 
magnetospheres. In the case of Jupiter, absorption by Callisto at ^ 27 Rj 
would be modified significantly by our results. However, no absorption 
signatures of this moon have been reported from either the Pioneer or Voyager 
particle experiments. At the orbit of Ganymede near 15 R,, where absorption 
has been seen [Burlaga et al,, 1980] along with signatures of a radially 
distended magnetic field CConnerney et al., 1981a], we have found the 
non-dipolar correction to be as large as a factor of 2.5 for equatorial 
particles and it should be incorporated into a quantitatively accurate 
analysis. For moons and rings of Jupiter lying within 14 R, of the planetary 
surface, effects of the magnetodisc lie mostly in modification of the bounce 
period (cf. Fig. 2a) for equatorial and near equatorial particles. 

In the case of Saturn, Rhea at 3.3 Rg lies in a region where significant 
non-dipolar effects occur. The deviation of Rhea's absorption as observed in 
the Voyager 1 cosmic ray experiment [Vogt et al., 1981] from its expected 
dipole position has been explained successfully using the ring current model 
CConnerney et al., 1981b], This suggests that the correction factors which we 
have calculated and which are most important in this radial region should be 
applied to any quantitative lifetime calculations. 

Finally we mention that because of the diminished scale length of the 
magnetic field near the equator, protons and heavier ions become non-ad iabatic 
at energies above a few hundred keV. This energy particle appears to dominate 
the particle energy density of the Jovian middle magnetosphere [Krimigis et 
al«, 1981] and may via its non-adiabatic meandering motion [Sonnerup, 1971] 
across the field reversal sheet contribute significantly to the magnetodisc 
current. 
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paraoittttrs of the Voyager MlddXe Magnetosphere Magnetic Models 


Jupite r Saturn 

*0 ^ »209 0 Rg^ 

Mglo 4,5X10"^ G 5X10"^ G 

Rq 5 Rj 8.5 Rg 

50 Rj 15.5 Rg 

2.5 Rj 2.5 Rg 
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FIGURE CAPTIONS 


Figs. 1. Hagnetio fiald lines in the model Jovian (Pig. la) and Saturnian 
(Fig. 1b) magnetospheres intersecting the equator at several values p^. 
Shown in dotted lines are dipole field lines which intersect the 
planetary surface at the same latitudes as model field lines (dipole p^ 
labels indicate the correspondence) . Note the different scales on the z 
and p axes. Because of strong magnetospheric plasma currents, Jovian 
field lines can be grossly non-dipolar; much of the Jovian middle 
magnetosphere maps to a small range of latitudes at the surface of the 
planet. Lines of constant magnetic latitude are also shown. 

Fi(^s, 2 , plots of the bounce period quantity H (cf. Eqs, 21 and 22) for 
particles in the model Jovian (Fig. 2a) and Saturnian (Fig. 2b) 
magnetospheres, H depends on the equatorial crossing distance of a 
field line and the particle's mirror latitude 1^^^ on that field line. As 
Pq 0, H approaches its value in a dipole field, which depends only on 
Ajjj. For each 1^^, the ratio HCp^, ijjj)/HC0 ,Xjjj) is the ratio of the bounce 
time of a particle mirroring at on a model field line passing through 
the equator at distance to the bounce time of a particle of the same 
energy mirroring at the same spatial point (p^, z^) but in the dipole 
field alone. 

Figs, 3. Plots of the azimuthal drift rate quantity F/G (cf. Eqs. 27 and 23) 
for particles in the model Jovian (Fig. 3a) and Saturnian (Fig. 3b) 
magnetic fields. F/G depends on the equatorial cro%^i^ing distance p^ of a 
field line and the particle's mirror latitude on that field line. As 
Pq 0, F/G approaches its value in a dipole field, which depends only on 

V’ ^m’ ^^^0' ^m^ is the 

ratio of the longitudinal drift rate of a particle mirroring at x^^^ on a 
model field line passing through the equator at distance p^ to the 
longitudinal drift rate of a particle of the same energy mirroring at the 
same spatial point (p^, z^) but in the dipole field alone. 

Figs. 4. Shown for the model Jovian (Fig. 4a) and Saturnian (Fig. 4b) 

-4 2 

magnetospheres is the quantity =3X10 B ( ZsO)/aB^/3z( z=0) which is 
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characttriatic In dtttroiinlng Units of adiabatlclty (of, Eq* 30). Shown 
also ar« values of proton kinetic energy as determined from Eq. (30) at 
which the proton gyro-radius equals the scale length of the field as 
determined from BizsQ)/iB^/iz(zsO) for representative values of 
equatorial pitch angle 6.. For electrons K.. is always greater tnan K. 
and hence the adiabaticity condition is less severe on them, as would be 
expected on the basis of their smaller rj^ss. 

Figs, 5. Shown for the model Jovian (Fig. 5a) and Saturnian (Fig, 5b) 

% 

magnetc^pheres is the value of electron kinetic energy at which the 
longitudinal drift rate <w> equals the planetary angular rotation rate 
^j,S* limited exceptions, which are discussed in the text, le 

sen.*.^ of <w> is prograde for protons and retrograde for electrons. 
Discontinuities in the Saturnian curves at 3.5 Rr. *>nd 15,5 R_ for 

w w 

equatorial particles reflect discontinuities in <u> owing to 

discontinuities in 3B_/3p which occur at the boundaries of the model 

z 

current sheet. At ./« 30 R , and p^ ^ Ro» For; equatorial 

particles because of the vanishingly small values of SB-/3P in these 
regions. Values of are approximately 1/2 to 3/4 as large as those 
shown here for electrons. 
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